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ABSTRACT 

The Gamma-ray Burst Monitor (GBM) on board Fermi allows to study the spectra of Gamma Ray Bursts (GRBs) over an unprece- 
dented wide energy range (8 keV - 35 MeV). We compare the spectral properties of short and long GRBs detected by the GBM (up 
to March 2010) with those of GRBs detected by the BATSE instrument on boai'd the CGRO. GBM and BATSE long bursts have 
similar distributions of fluence (F), E'^^^^ and peak flux (P) but GBM bursts have a slightly harder low-energy spectral index a with 
respect to BATSE GRBs. GBM and BATSE short bursts have similar distributions of fluence, a and peak flux, with GBM bursts 
having slightly larger E'^^^^. We discuss these properties in light of the found correlations between -Ep^^ fluence and the peak 

flux. GBM bursts confirm that these correlations are not determined by instrumental selection effects. Indeed, GBM bursts extend the 
E'^^-F and E'^^^^-P correlations both in fluence/peak flux and in peak energy. No GBM long burst with E°^^^_ exceeding a few MeV 
is found, despite the possibility of detecting it. Similarly to what found with BATSE, there are 3% of GBM long bursts (and almost 
all short ones) that are outliers at more than 3(T of the iipeak " £^iso correlation. Instead there is no outlier of the fip^ak - ^iso correlation, 
for both long and short GBM bursts. 

Key words. Gamma-ray burst: general - Radiation mechanisms: non-thermal 



1. Introduction 

The Fermi satellite, launched in June 2008, offers a great op- 
portunity to characterize GRB spectra over a wide energy range 
thanks to its two high-energy instruments: the Large Area 
Telescope (LAT) and the Gamma-ray Burst Monitor (GBM - 
Meegan et al. 2009). The GBM is composed by twelve Nal de- 
tectors (with good spectral resolution between ~8 keV and -1 
MeV) and two BGO detectors (operating between 200 keV and 
40 MeV). Significant emission in the LAT energy range (~ 30 
MeV - 300 GeV) has been detected only in about 20 GRBs 
till now (December 2010), while the GBM triggered about 600 
GRBs. 

A detailed spectral analysis of all GRBs detected by the 
Fermi/GBM up to the end of Mai'ch 2010 (438 events) has been 
performed by Nava et al. (2010; NIO hereafter). These spectra 
were fitted with different models and for 316 events (272 long 
and 44 short) it was possible to constrain the peak energy of the 
vFy spectrum (£°g^i^). For long bursts we found {E'^l^} ~ 160 

keV and an average low-energy power law index {a} 0.9. 

Short bursts are found to be harder, both in terms of £'°''^ and a: 

peak 

<£peak) ~ 490 keV and (a) ~ -0.5. 

NIO also analyzed the peak spectrum of GBM bursts, i.e. 
the spectrum corresponding to the peak flux of the light curve, 
accumulated on a timescale of 1.024 s and 0.064 s for long 
and short events, respectively. The comparison with the time- 
integrated spectral properties shows that peak spectra, on aver- 
age, have harder low-energy spectral indices but similar peak 
energies with respect to time-integrated spectra. 
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Before Fermi, other instruments allowed to study the prop- 
erties of the prompt emission spectra of GRBs. Due to its broad 
energy range (~25 keV 2 MeV), high sensitivity and detec- 
tion rate, the BATSE instrument onboard the Compton Gamma 
Ray Observatory ( CGRO) satellite has been so far the instrument 
best suited to characterize the GRB prompt emission properties. 
Thanks to its almost all-sky viewing, BATSE detected more 
than 2700 GRBs in about 9 years. 

The published spectral catalogs of BATSE bursts comprise 
relatively small sub-samples of bright GRBs, selected on the 
basis of the burst fluence and/or peak flux (Preece et al. 2000; 
Kaneko et al. 2006 - K06 hereafter). The analyzed samples al- 
lowed to study the spectral properties only of long bursts, given 
the small number of short bursts present in these samples (e.g. 
17 short GRBs in the K06 sample). These studies revealed that 
the low-energy power law index distribution of long GRBs is 

centered around a 1 and pointed out the inconsistency of 

the large majority of burst spectra with a synchrotron interpre- 
tation. The distribution of long GRBs analyzed by K06 

peaks around 250 keV, with a relatively narrow disper- 

sion. However, this refers to bright bursts, while fainter bursts 
have smaller E°^\ values, as shown by Nava et al. (2008, N08 

peak ' ^ ' 

thereafter). They performed the spectral analysis of a sample of 
BATSE bursts selected by extending the limiting fluence of K06 
(i.e. F - 2x 10"'' erg cm"^) down Xo F - 10"'' erg cm"^. They 
found that £'°''^ correlates with the fluence F and the peak flux 

peak ^ 

P. This sample of BATSE faint bursts has a distribution of £'°'"*, 

^ peak 

values centered at ~ 150 keV, i.e. a value smaller than the one 
found for the bright BATSE bursts analyzed by K06, as a con- 
sequence of the mentioned E°^\-F correlation. This result con- 

^ peak 
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firms that the derived distribution of E'^^^^ is strongly affected by 
the adopted cuts in fluence (or peak flux). 

N08 also found a correlation between E'^^^^ and the flu- 
ence/peak flux for short bursts. This implies that when we com- 
pare the £pg2]j distributions of short and long GRBs we must take 
into account the possible different fluence/peak flux selection 
criteria. A large sample of short BATSE bursts have been ana- 
lyzed by Ghirlanda et al. 2009 (G09 hereafter). They performed 
a detailed spectral analysis of 79 short bursts and compared their 
properties with those of 79 long BATSE bursts selected with the 
same limit on the peak flux. They found that the E°^\ distribu- 

^ peak 

tions of the two classes are similar, while the low-energy power 

law indices are different: short bursts have (a) 0.4, harder 

than long events. 

A well known property of GRBs, related to their prompt 
emission, is the correlation between the rest frame peak en- 
ergy fipeak of long bursts with the bolometric isotropic energy 
iiiso emitted during the prompt (Amati et al. 2002) or with the 
bolometric isotropic luminosity Lp^so estimated at the peak of 
the light curve (Yonetoku et al. 2004). Such correlations rep- 
resent an intriguing clue on the dominant emission mechanism 
of the prompt phase. Furthermore, if corrected for the jet open- 
ing angle, their dispersion reduces considerably (Ghirlanda et al. 
2004a) and allows to use GRBs as standard candles (Ghirlanda 
et al. 2004b). 

The correlations in the observer frames (^p^^-^ and ^pg^i^- 
P) may be just the consequence of the rest frame (iSpeak - £^i,so 
and /ipeak - Lisa respectively) correlations mentioned above. 
Alternatively, it has been claimed that the rest frame correlations 
are the result of instrumental selection effects (Band & Preece 
2005; Nakar & Piran 2005). Ghirlanda et al. (2008, hereafter 
GOB) and N08 examined the instrumental selection effects which 
may affect the observer frame correlations E°^\ -F and E°^'^, -P. 

peak peak 

They found that, although instrumental biases do affect the burst 
sample properties, they are not responsible for the correlations 
found in the observational planes. 

Moreover, Ghirlanda et al. (2010a) recently showed that 
the correlation fipeak-^iso and iSpeak-^'p.iso holds for the time- 
resolved quantities within individual long GBM bursts (see also 
Firmani et al. 2009 for Swift bursts and Krimm et al. 2009 for 
Swift-Suzaku GRBs) and that this "time-resolved" correlation is 
similar to that defined by the time-integrated properties of dif- 
ferent GRBs. Similar results were found for short GBM bursts: 
there is a significant correlation between the observer frame peak 
energy E'^^^ and the peak flux within individual short GRBs 
and this correlation has a slope similar to that of the rest frame 
£^peak - iiso Correlation (Ghirlanda et al. 2010b). These results 
confirm that the "Amati" and "Yonetoku" correlations have a 
physical origin, instead of being the result of instrumental selec- 
tion biases as claimed, and that the trends {E°'°^, -F and E°'°'^, -P) 

' ^ peak peak ^ 

seen in the "observational planes" are just their outcome. 

Through the spectral catalog of GBM bursts of NIO, we can 
study the distribution of GBM bursts in the observational planes 
E°\-F and E°\-P and test if the coiTelations found by BATSE 

peak peak 

bursts (G08, N08 and G09) still hold. In the observational planes 
we can also study, for the first time, the possible instrumental 
biases of GBM for the bursts analyzed in NIO, and compare the 
spectral properties of long and short GRBs detected by BATSE 
and by the GBM. Finally, we can also compute the fraction of 
GBM short and long GRBs which are outliers, for any assigned 
redshift, of the rest frame fipeak-^^iso and iipeak-i'iso correlations. 
These are the main aims of the present paper. 



In §2 we present the samples of BATSE and GBM bursts 
(both long and short) used for our comparison. We compute 
(§3) the relevant instrumental selection effects introduced by the 
GBM on the observational E'^^^-P and E'^^^-F planes consid- 
ering separately short and long GRBs. The comparison between 
BATSE and GBM results is also presented in terms of spectral 
parameters distributions in §4. We discuss our results and draw 
our conclusions in §5. 

2. Samples 

2.1. Long bursts 

BATSE — Fig. [Ushows the LogA^-Logf of long GRBs detected 
by BATSE (open squares), where is the number of objects 
with fluence larger than F. Since we compare this LogA^-LogF 
with the one of GBM bursts (filled squares in Fig. [TJ, we com- 
pute F in the energy range between 20 keV and 2000 keV, which 
is common to both instruments. BATSE fluences are taken from 
the online CGRO/BATSE Gamma-ray Bursts Catalo£l 

We have then considered the spectral catalog of K06 which 
contains aU BATSE bursts with peak flux P (50 - 300 keV) > 
10 photons cm-^ s ' or fluence F (20 - 2000 keV) > 2 x 
10~^ ergs cm"^. From the 350 events of this sample we extract 
the long ones, i.e. with observed duration Tgo > 2 s. Among 
these, 280 GRBs have a well determined (i.e. their spec- 
tra are best fitted by a curved model with a < -2 or /? > -2, 
i.e. showing a peak in a vFy representation). For 104 events the 
spectrum is best fitted by a BAND model (Band et al. 1993), for 
65 by a COMP model (a power law with a high-energy exponen- 
tial cutoff) and for 1 1 1 by a smoothly broken power law model 
(SBPL). 

Although the K06 analysis enlarges the previous spectral cat- 
alog of BATSE bursts (Preece et al. 2000) it still selects only the 
brightest BATSE bursts (i.e. corresponding to only the 13% of 
the entire population of bursts detected by BATSE). N08 se- 
lected a sample of 100 BATSE bursts with a fluence fainter 
than the threshold adopted by K06. The N08 events, in fact, 
have a fluence in the range 10"^ < f < 2 X 10"^ erg cm . 
Moreover, although these are only 100 bursts, they are repre- 
sentative of the large population of ~1000 bursts in this fluence 
range since they were randomly extracted following the LogA^- 
LogF distribution of BATSE GRBs in this fluence range. Of 
these 100 representative bursts, 44 are best fitted by the COMP 
model, 44 by the BAND model, and 12 by a power law (PL) 
function. Therefore the N08 sample contains 88 GRBs with a 
spectrum fitted by a curved model, for which £'°''^ was well de- 

^ ' peak 

termined. Among the bursts analyzed by K06 with peak flux P 
larger than 10 ph cm"^ s ' , there are GRBs with fluences smaller 
than F - 2y. 10"^ erg cm""^, that overlap with the ones studied 
by N08. We exclude these bursts from the present discussion, 
in order to have well defined complete samples at two limiting 
fluences, that we will call, in the rest of the paper, the "bright" 
BATSE bursts (the bursts in K06 with F > 2y.\Q-^ erg cm"^) 
and the "faint" BATSE bursts (the bursts studied in N08 with 
10"^ < F < 2x 10"^ erg cm"^). The grey shaded regions in Fig. 
[T]correspond to this subdivision. 

GBM — In NIO we have analyzed the spectra of all the GRBs 
detected by the GBM up to March 2010 (438 GRBs). No fluence 
or peak flux selection has been adopted. Fig. [T] shows that the 
shape of the LogA^-LogF for the two instruments is very similar. 
To compare the E°^^^ and a distribution of GBM bursts with 
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Fig. 1. Long GRBs. LogA^-LogF for long BATSE bursts (empty 
squares) and long GBM bursts (filled circles). In both cases the 
fluence F is integrated between 20 keV and 2000 keV. GBM 
data are taken form NIO, while BATSE data are from the online 
catalog and include all the BATSE bursts for which the fluence 
has been estimated. For reference two power laws with slope - 
3/2 are shown (dot-dashed lines). 



those of BATSE bursts, we select from the NIO catalog two 
subsamples with the same fluence criterion adopted by K06 and 
N08 for BATSE bursts and we call them the bright GBM sample 
and the faint GBM sample. 

In six GBM bursts NIO could not analyze the spectrum, due 
to lack of data. For the remaining 432 bursts we performed the 
spectral analysis using a power law model (PL - 110 spectra), a 
COMP (232 spectra) and a Band model (90 spectra) and evalu- 
ating for each burst the spectral parameters of the best fit model. 
359 events belong to the long burst class. We also estimated their 
peak flux on time bin of 1.024 seconds. In this work we will use 
this sample of GBM long bursts for comparison with BATSE 
long GRBs. 

2.2. Short bursts 

BATSE — The most comprehensive sample of short BATSE 
GRBs with well defined spectral parameters is composed by 
the 79 events analyzed by G09, selected for having P > 
3 photons cm"^ s"'. In 71 cases the spectra have a well deter- 
mined E'*'^ 



^peak' 



GBM — For the GBM instrument we use the spectral pa- 
rameters of the 44 short GRBs present among the 438 bursts 
analyzed by NIO with a well determined ^^^i^. Their peak flux 
is estimated on a time bin of 0.064 seconds. 



3. £°''\-Fluence and £°''\-Peak Flux planes: 

peak peak 

comparison between BATSE and GBM bursts 

A correlation between the total fluence and E'^^^^ was first found 
by Lloyd, Petrosian & Mallozzi (2000) for a sample of BATSE 
bursts without measured redshifts. This finding was recently 
confirmed by Sakamoto et al. (2008) using a sample of bursts 
detected by Swift, BATSE and Hete-11. In particular, they noted 



that X-Ray Flashes and X-Ray Rich bursts satisfy and extend 
this correlation to lower fluences. 

The distribution of GRBs with and without measured red- 
shift in the planes E'^^^^-F and E'^^^^-P has been investigated by 
G08 and N08. N08 considered all events with published spectral 
information detected by different instruments (Swift, BATSE, 
Hete-II, Konus/Wind and BeppoSAX) together with the 100 
faint BATSE bursts analyzed in that paper. In both planes long 
bursts define a correlation, with fainter bursts having lower E'^'^^. 

In order to examine the distribution of GBM bursts in the ob- 
servational planes E°^'^, -F and E°^'^, -P and compare it with the 

^ peak peak ^ 

BATSE bursts we have first to estimate the possible instrumen- 
tal biases induced by the detector (see G08). One instrumental 
bias is the capability of an instrument to be triggered by a burst, 
i.e. the "trigger threshold" (TT) (first computed by Band 2003 
for different detectors). The second bias concerns the minimum 
number of photons required to analyze the spectrum and con- 
strain the spectral parameters. This is called "spectral threshold" 
(ST) in G08 and N08. The TT translates into a minimum peak 
flux, which depends on the burst spectrum and in particular on its 
E°^'', and can be described as a curve in the E°^\ -P plane. The 

peak peak ^ 

second requirement (ST) results into a minimum fluence which 
depends on E^^ and also on the burst duration. For this reason 

the ST is represented as a region (i.e. not a line) in the E'^^^^-F 
plane. 

These curves (TT and ST) divide the observational planes 
^peik~^ and E'^^^-P into two regions. Bursts with peak energy 
and peak flux which place them on the left of the TT curve 
cannot be triggered by the corresponding instrument. Similarly, 
bursts with peak energy and fluence which place them on the 
left side of the ST curves do not have enough photons to allow a 
reliable spectral analysis (see G08 for more details). 



3.1. Estimate of GBM instrumental selection effects 

Following G08, the TT curves are obtained adapting the results 
of Band (2006) and are shown in the right (top and bottom) pan- 
els of Fig. |2] The TT curves are the same for long and short 
bursts as the trigger threshold depends only on the peak flux. 

The ST curves have been calculated from numerical simu- 
lations, as described in G08. To perform these simulations, the 
typical background spectrum and the detector response function 
are required. For the bursts detected by the GBM both of them 
depend on several factors (e.g. the satellite attitude when a burst 
occurs), and thus there is no universal background and response 
matrix which can be adopted. To overcome this problem we use 
the real backgrounds and responses of several GRBs detected 
by the GBM and average the results of the simulations to build 
average ST curves for the population of long and short bursts, 
respectively. 

Our simulation performs a joint spectral analysis of spectra 
simulated for two Nal detectors and one BGO detector For long 
bursts, simulations were performed using the detector response 
files and the background spectra of the long GBM bursts pub- 
lished in GIO. We considered, as done in G08, two representa- 
tive values of the duration i.e. Tgo = 5 and 20 s (corresponding 
in Fig. |2]to the curves delimiting the red shaded region on the 
left and right side, respectively). For short bursts we estimated 
the ST curves adopting the response files and the background 
spectra of the short bursts of the GBM sample of NIO assuming 
a typical duration of the simulated spectra of 0.7 s (red curve in 
the bottom left panel of Fig. [2]). This value, also adopted by G08 
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Fig. 2. /ip^^i^-Fluence and /s^^j^-Peak Flux planes for long (upper panels) and short (bottom panels) bursts. Empty squares represent 
BATSE bursts, filled circles GBM bursts and filled triangles indicate events detected by other instruments (from N08). In all panels 
the instrumental limits for BATSE and GMB are reported: shaded curved regions in the upper left panel show the ST, estimated 
assuming a burst duration of 5 and 20 s; solid curves in the bottom left panel represents the ST for short bursts. Solid curves in the 
right panels define the TT, identical for short and long events. Thresholds for BATSE are taken from G08 while those for the GBM 
instrument are derived in this work (see §3.1). The dashed curve in the bottom right panel represents the selection criterion applied 
by G09 for their sample of short bursts, i.e. P > 3 phot cm"^ s The shaded regions in the upper left corners of all the planes 
are the region identifying the outliers at more than 3cr of the Zipeak - £^iso (left panels) and Zipeak - ^^iso (right panels) correlations 
for any given redshift. GRBs, without measured redshift, which fall in these regions are outliers of the corresponding rest-frame 
correlations (fipeak - £^iso and iipeak - hso for the left and right panels respectively) for any assigned redshift. It means that there is 
no redshift which make them consistent with these correlations (considering their 3cr scatter). 



for BATSE bursts, corresponds to the typical duration of short 
GRBs observed by the GBM. 

For the TT and ST curves of the BATSE instrument we sim- 
ply report those obtained in G08 (for long bursts) and in G09 
(for short bursts). 

Fig. |2] shows the distribution of GBM bursts in the ^p^^i^- 
F and E'^^^^-P planes (right and left panels) for long and short 
GRBs (upper and bottom panels). 



3.1.1. Long bursts 

^peak Fluence — As can be seen in the top left panel of Fig.|2] 
the distribution of long GBM bursts (filled circles) extends down 
to the lower end of the distribution of BATSE bursts (empty 
squares). 

The presence of GBM bursts with low E'^^^^ (between ~10 
and ~50 keV), not present in the BATSE sample, is clearly due to 
the wider energy range of the GBM instrument, sensitive down 
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to ~8 keV (see the ST curves). In this region, GBM bursts are 
consistent with bursts detected by other instruments (filled trian- 
gles). We also note that GBM bursts define a correlation which 
mostly overlaps with that defined by BATSE bursts and extends 
to the lower-left part of the E'^^-^-F plane. 

Despite the GBM assures good coverage up to ~30 MeV (vs 
~ 1 MeV of BATSE), there are only a few long GRBs with E'^l^. 
exceeding few MeV, similarly to what found for the population 
of BATSE bursts. Note that high E°^l^. also means GRBs with 
high fluences, which are rarer than GRBs with low fluences. 

Note also that while BATSE bursts appear concentrated at 
the high end of the E'^^^^^-F correlation, the sample is composed 

by all bursts with F > 2x 10"^ erg cm"^ (the K06 sample) and 
only one hundred of fainter bursts (the N08 sample) represen- 
tative of ~ 1000 objects with 10"'' < F < 2 x 10"^ erg cm"^. 
Therefore the real density of BATSE bursts in the latter fluence 
range is much larger than what represented in Fig. |2] (see e.g. 
Fig. 8 in N08) so that the sHght shift between BATSE and GBM 
population density in the upper panels is only an apparent effect. 

Another result shown by the GBM bursts and consistent with 
the conclusion drawn from BATSE is that the ST effect is not 
responsible for the distribution of the data in the plane, i.e. it 
cannot explain why bursts tend to distribute along a correlation. 
This is well visible for BATSE bursts: events with large E°^\ 

^ peak 

tend to concentrate far from the ST, i.e., at higher fluences. The 
trend shown by the BATSE ST (which requires higher limiting 
fluences when ^^^i^ is very high and very low) cannot explain 
this behavior The same holds for GBM bursts, for which the ST 
are even less curved and cannot be responsible for the observed 
correlation which has a slope ~ 0.20 + 0.04, consistent with that 
defined by BATSE bursts (N08). 

^peak P^<^k Flux — The distribution of GBM long bursts in the 
E'^l^-P plane with respect to BATSE bursts is shown in the up- 
per right panel of Fig. |2l Solid curves represent the TT derived 
for both instruments (adapted from Band et al. 2006). On aver- 
age, the GBM instrument is a factor 3 less sensitive than BATSE 
in the common energy range. As remarked by N08, the sample 
of BATSE bursts lies far from its TT, suggesting that for this in- 
strument the demand of performing a reliable spectral analysis is 
the dominant selection query. This is not the case for the GBM: 
the data points lie very near the TT curves, suggesting that if a 
burst is detected there is a good chance to recover its spectral 
parameters. This makes TT and ST competitive selection effects 
for GBM bursts. 



3.1.2. Short bursts 

For short bursts (bottom panels in Fig. |2]) the situation is quite 
different in both E°^'', -F and E°^\ -P planes than for long 

peak peak ^ ^ 

events. 

Although GBM short bursts are still only a few (filled cir- 
cles in Fig. |2]) there is a weak indication of a correlation in 
both planes (left and right bottom panels in Fig. |2l), consistent 
with the trends suggested by BATSE short bursts (open squares). 
However, the overall behavior is different from what happens for 
long GRBs. 

At high GBM short events occupy the same region of 
the BATSE ones, and even extend the E°^\-F trend to E°\ val- 

' peak peak 

ues larger than 1 MeV (i.e. above the BATSE upper threshold). 



revealing that short GRBs with E^^^^ larger than 1 MeV exist in 
the population of GBM events. 

Furthermore, contrary to what one might expect, there are no 
short GBM bursts wifli E°^l^ below ~ 200 keV. At low fluences 
this can be accounted for by considering the ST derived for the 
GBM instrument: Fig. |2] shows that GBM short events lie very 
near to the ST curve (bottom left panel of Fig.|2]i, that prevents 
the estimate of E°ll^,< 200 keV when F < 4 x lO^^erg cm"-. 

The higher sensitivity of BATSE instead implies that its ST is 
located at lower fluences. This however does not account for the 
absence of short GBM burst with fi*', < 200 keV and fluence 

peak 

> 5 X 10"^ erg cm"^: the BATSE sample shows that GRBs with 
low £'°''', but on the right side of the GBM ST do exist. Their 

peak ^ 

absence in the present GBM sample seems to imply that they 
are relatively rare and therefore that at a fluence > 4 x 10"^ erg 
cm"^ most of the events should have £pg^|j> 200 keV. In other 

words this would support the presence of a E'^^^-F correlation 
for short GRBs albeit with a large dispersion. A larger sample of 
short GBM bursts is required to confirm this. 

In terms of E°^\ this translates into a distribution peaked 

peak ^ 

at higher energies compared both to the BATSE one and to the 
GBM long events (see §. 4.2). 

In the ^peaij-^* plane of short GRBs (right bottom panel in 
Fig. |2| similar conclusions can be drawn: for both instruments 
the TT curves (solid lines) do not affect the samples. Both sam- 
ples are clearly limited by the selection cut applied on the peak 
flux (shaded curve). From the fact that the peak flux of GBM 
bursts is significantly above their TT, we infer that their selec- 
tion is dominated by the ST. 



3.2. Outliers of the Ep^ak - -Eiso and ispeak - ^iso correlations 

Another relevant point is to test whether bursts without measured 
redshifts are consistent with the fipeak - £^iso and fipeak - ^iso cor- 
relations. These correlations are defined in the rest frame and re- 
quire, to add a burst on top of them, to have the redshift known. 
However, as first proposed by Nakar & Piran (2005) and then by 
Band & Preece (2005), knowing ^p^^k ^^'^ fluence or peak 
flux it is still possible to test if a burst, without measured red- 
shift, is an outlier of the ZSpeak - ^iso and fipeak - i^iso correlations. 

In the observational planes it is possible to define a "region 
of outliers". We start by writing the fipeak-^^iso correlation as (the 
same argument can be repeated for the fipeak - Li^o correlation): 



Epeak - KE^^^ 



(1) 



Since E°°l^^^ £peak/(l + z) and F = E^Ul + z)/4;rc/^(z), we can 
form the ratio. 



' peak^ 



Andliz) 



(l+z)( 



1 +'?)/'? 



(2) 



where cr corresponds to the scatter of data points around the rest 
frame correlation that is being tested. Note that this is not the 
error on the slope or normalization of the correlation, but the 
scatter measured perpendicular to the best fit line of the iipeak 
/ correlation and modeled as a gaussian distribution. 

The RHS of the above relation is a function of z, rj and 
cr only. The upper limit of the ratio E^^J'^^^^^/F establishes an 
allowance region boundary on the corresponding plane (upper 
left-corner shaded regions in Fig. |2]i. All the bursts that fall 
below this line in the observational planes can have a redshift 
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which make them consistent with the iSpeak - ^iso coirelation 
within its 3cr scatter. Those falhng in the shaded region are out- 
liers at more than 3cr for any assigned redshift. 

Ahhough this test has been akeady appHed several times in 
the recent past, some guidelines should be followed: 

1. since the rest frame correlations are defined with the bolo- 
metric 1 keV-10 MeV Eiso and Lp^so, then when testing the 
region of outliers in the observational planes one should use 
the fluence F and peak flux P defined on the same energy 
range; 

2. it is correct to consider the 3cr scatter of the rest frame corre- 
lations and not the uncertainty on the slope {tj) and normal- 
ization (K) of the correlations. This is because the scatter cr 
of the Zipeak-^^iso and iipeak-i'iso Correlations dominates over 
the statistical uncertainty on K and 77 (e.g. GIO); 

3. while the fipeak - £^iso and /Speak - Lino correlations were first 
derived with only a dozen of bursts, they have been now up- 
dated with nearly 100 GRBs with measured redshifts (e.g. 
N08, GIO): the correlation parameters (slope, normaliza- 
tion and scatter) have changed since their discovery, so one 
should adopt the most updated versions of these correlations. 

N08 find that 6% of BATSE long bursts are outliers of the 
Epe-ik - £iso correlation, while no outlier is found for the fipeak- 
Lp;Ko correlation. Almost all short BATSE bursts are outliers of 
the /ipeak-^^iso Correlation (defined by long bursts) but they can 
be consistent with their very same /Speak-i-p.iso correlation. These 
findings are supported by the consistency of the few short bursts 
with measured redshift with the fipeak - hso correlation while 
they are outliers at more than 3cr of the fipeak - ^iso correlation 
(G09). 

For GBM long bursts. Fig. |2] shows that only 8 GRBs (i.e. 
3%) lie in the region of outliers (in the ^p^^^-^^ plane), even if 
the extended energy range of the GBM allows to explore the 
region of large E'^^^^^ and intermediate fluences, where outliers, 
if they exist, could be found. 

Instead, for short GBM bursts, we can see that most of them 
are outliers at more than 3cr of the iSpeak - £^iso correlation. On 
the contrary, there are no short bursts in the region of outliers 
in the E'^^-^-P plane, i.e. they are all consistent with the fipeak - 
Liso correlation defined by long GRBs. The hypothesis that short 
and long bursts follow the same iSpeak-ip.iso correlation is also 
supported by the few short events with known redshift (G09). 



4. Spectral parameter distributions 

In this section we compare the distributions of the spectral pa- 
rameters (low energy spectral index a and peak energy ^pg^j^) 
for long (§4.1) and short (§4.2) GRBs detected by BATSE and 
by the GBM. In addition, we show the LogA^-LogP distribution 
for short bursts. 



4.1. Long Bursts 

Fig. |3] (upper panel) shows the distributions of E'^^^^^ of bright 
BATSE and bright GBM bursts. The two distributions are quite 
similar: the central value of the gaussian fit (solid line) is 
E°\~260 keV and £'*V~280 keV, respectively, with GBM 

peak peak ^ r j ^ 

bursts having a larger distribution (standard deviation cr=0.33 
to be compared with cr=0.21 for BATSE bursts) and extending 
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Fig. 3. Long GRBs. E°^J^^ distribution for BATSE bursts (blue 
solid filled histograms) and GBM bursts (red hatched his- 
tograms). For both instruments, we plot separately the £pg^k 
tributions for the bright sample (upper panel) and the faint sam- 
ple (bottom panel). 



Using the Kolmogorov-Smirnov (KS) test, we find that the prob- 
ability that the two distributions are drawn from the same parent 
population is 0. 164. 

For faint bursts (bottom panel of Fig. [3]) the results are simi- 
lar. For both instruments the gaussian fit to the E°^^^^^ distribution 
is centered around 140 keV and has a standard deviation cr=0.31 
(the KS test probability is 0.18). Also in this case the GBM dis- 
tribution is larger In particular, from the comparison of the two 
histograms it is appears that GBM data allows to recover very 
low £pg^^, thanks to the good GBM/Nal sensitivity down to 8 
keV (i.e. extending by a factor of 3 the low-energy bound of 
BATSE). For BATSE bursts there is a quite sharp cutoff' at ~ 
50 keV. This is in agreement with the simulations performed by 
G08 (see Fig. |2]i showing that it is very difficult for BATSE to 
recover E'^^^^ < 50 keV. Large fluences, unusual for such values 

of £pg3|,, would be required. All the results of the gaussian fits 
and KS probabilities are summarized in Tab. [T] 

We stress that the E'^^^^^ properties of faint and bright bursts 

are very different due to the correlation between E°'°.l,, and F: the 



both at lower and higher /i^^j^with respect to that of BATSE. 



faint sample is characterized by a central value of £'°''!,, which is 
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Fig. 4. Long GRBs. Low-energy power law index a distribu- 
tions for BATSE bursts (blue solid filled histograms) and GBM 
bursts (red hatched histogram). For both instruments, we plot 
separately the distributions for the bright sample (upper panel) 
and the faint sample (bottom panel). For the bright BATSE sam- 
ple of K06 (upper panel) we also show separately the a distri- 
bution of bursts fitted with all models (light gray shaded his- 
togram) and that of bursts fitted with only the COMP or Band 
model (dark blue histogram). 



almost a factor of 2 lower than that of the bright sample. The KS 
probability of the distributions of E'^^y^ within the GBM sample 

between faint and bright bursts is 5.7x10"^. 

Fig. |4] shows the a distribution for bright (upper panel) and 
faint burst (bottom panel). In both cases, GBM bursts tend to 
have a harder low-energy power law index and a somewhat 
tighter distribution (central values and standard deviations of the 
gaussian fits are reported in Tab.[TJ. However the KS test shows 
that the a distributions of faint GBM and faint BATSE bursts are 
similar (KS probability=0.12). Moreover, GBM bursts do not 
show a significant relation between a and the fluence because 
faint and bright GBM bursts have similar distributions peaked 
respectively at -0.91 and -0.94 and with a KS probability of 
0.03. 

The a distribution of bright BATSE bursts (grey histogram), 
instead, shows significant difference with respect to both the 
bright GBM sample (KS probability=6xl0"^) and the faint 
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Fig. 5. Short GRBs. LogA^-LogP for short BATSE bursts 
(empty squares) and short GBM bursts (filled circles). In both 
cases the peak flux P (in photons cm"^ s"') is integrated be- 
tween 50 keV and 300 keV on a time scale of 64 ms. GBM data 
are taken form NIO, while BATSE data are from the online cata- 
log and include all the BATSE bursts for which the peak flux has 
been estimated. The vertical dashed line is the flux limit of the 
selection of BATSE bursts analyzed in G09 and for reference a 
power law with slope -3/2 is shown (dot-dashed line). 



BATSE sample (KS probability=7xl0"^). Bright BATSE bursts 
tend to have softer a with respect to all the other samples. We in- 
vestigated the possible origin of this difference considering that 
the K06 sample contains bursts whose spectra are fitted with the 
SBPL, COMP or Band model. As noted by K06 themselves the 
spectral parameters E'^^-^ and a do show different typical values 
and width of their distributions depending on the fitting spectral 
model. Considering that GBM bursts are adequately fitted by ei- 
ther the COMP or the Band model, we excluded from the ^p^^i^ 
distribution of BATSE bursts the bursts fitted with a SBPL. The 
resulting histogram (dark blue shaded in Fig.|4|i is now fully con- 
sistent with the distribution of a of GBM bursts (the KS proba- 
bility now becomes 0.4). 



4.2. Short Bursts 

The largest sample of BATSE bursts for which the spectral anal- 
ysis has been performed was selected on the basis of a peak flux 
criterion (G09). A meaningful comparison between GBM and 
BATSE short bursts requires a sample of short GBM bursts se- 
lected on the basis of the very same criterion. Before investigat- 
ing the spectral parameter distributions, we compare the LogA^- 
LogP for both instruments, where P in this case is the peak flux 
in photons cm"^ s Since for the GBM sample NIO estimate 
the peak flux on the 64 ms timescale, also for BATSE bursts we 
select (from the online catalog all bursts for which the P on 
64 ms has been estimated. For BATSE bursts the peak flux is 
integrated in the 50 keV - 300 keV energy range. Therefore, we 
estimate for all short GBM bursts in NIO the photon peak flux 
between 50 keV and 300 keV. Fig.|5]shows our results. 

As discussed in §3 the lack of bursts with low peak flux in 
the GBM sample is due to the ST threshold shown in Fig.|2] This 



http://heasarc.gsfc.nasa.gov/W3Browse/cgro/batsegrb.html 
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Fig. 6. Short GRBs. fi™^^ distribution for BATSE bursts (green 
hatched histogram, from G09), and GBM bursts (purple filled 
histogram, from NIO). 
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Fig. 7. Short GRBs. a distribution for BATSE bursts (green 
hatched histogram, from G09), and GBM bursts (purple filled 
histogram, from NIO). 



instrumental threshold, indeed, dominates over the TT thresh- 
old and determines that short GRBs for which the spectrum can 
be analyzed and the spectral parameters properly constrained 
should have a large number of photons. At high peak fluxes in- 
stead, GBM has detected more short GRBs than BATSE due 
to the ^pggij-^' correlation, which associates high peak fluxes to 
high peak energies, the latter better constrained with the larger 
energy range of the GBM instrument (BGO) than with BATSE. 
This explains the different shapes of the two LogA^-LogP. 

To compare the spectral parameters, we consider short GRBs 
(from G09 for BATSE and NIO for GBM) with known a and 
and with P(50-300 keV)>3 ph cm^^ The dis- 
tributions are shown in Fig. |6] and are quite different. The lack 
of low in the GBM sample (which corresponds to the lack 
of low peak fluxes in Fig. |6l) can be explained by considering 
the shape and position of the ST for short bursts (left bottom 
panel in Fig.|2]l and the existence of a correlation between P and 
^peak ^^^^ discussion above and in §2). Moreover, contrary to 
long bursts, the E°^^^^ distribution of GBM short events extends 
to higher energies, suggesting that such large values of ^p^^i^ can 
be found in short bursts and that they were not present in the 
BATSE catalog due to its limited energy range (up to only ~ 1 
MeV). 

The a distributions of BATSE and GBM bursts is con- 
siderably different. The GBM confirms that short events are 
harder than long ones in terms of low-energy spectral index 
(KS probability=3xlO"^). However, the distribution is peaked 
around a - -0.59 and it is very narrow (cr=0.15). We tentatively 
interpret this as due to the large energy range of GBM which ex- 
tends down to 8 keV, but this point deserves further study. The 
extension down to low energies of GBM allows in principle to 
determine a more accurately. Instead, the limited energy range 
of BATSE resulted in less accurate estimates of a and thus a 
more dispersed distribution of its values. This effect is slightly 
present also in long bursts, both bright and faint (see Tab.[TJ. 



5. Discussion and Conclusions 

In this work we investigated the presence of the £'°^^j,-Fluence 
and /ip^^j.-Peak Flux correlations in GBM bursts (both long 
and short) detected by the GBM instrument up to March 2010. 
Similarly to what has been done for long and short GRBs de- 
tected by BATSE (N08; G09) we examined the distribution of 
GBM bursts in the £'°''^, -Fluence and E°^\ -Peak Flux planes in 

peak peak ^ 

order to study instrumental selection effects and test their con- 
sistency with the rest frame correlations (i.e. fipeak - £^iso and 
£^peak - Lisa, respectively) defined by GRBs with measured red- 
shifts. To this aim, we have estimated, for the GBM instrument, 
the spectral threshold (ST) and the trigger threshold (TT), in 
order to quantify the selection effects acting on the considered 
samples and their role on the found correlations. 
Our main results are: 

- Long GRBs detected by GBM follow the same E'^l^_-F and 
^peak"^ correlations defined by BATSE GRBs (Fig.|2]i. We 
computed the instrumental selection effects of GBM - as al- 
ready done for BATSE (G08; N08): the trigger threshold and 
the spectral threshold are not responsible for the correlations 
defined by long GRBs in both planes (see Fig.|2]i. The GBM 
spectral extension down to 8 keV with respect to the limit of 
30 keV of BATSE, allows to extend the correlations to lower 
peak energies/fluences. Instead, despite of the higher energy 
threshold of GBM (40 MeV) no long GRB with E°ll^. larger 
than a few MeV is detected. This can be due to a real ab- 
sence of bursts with such high ^p^^i^ or to the fact that they 
have large fluences, thus being too rare to be detected during 
less then 2 years of the GBM observations. 
We conclude that long GRBs detected by GBM confirm what 
found with BATSE bursts, in particular that they follow a 
correlation both in the E'^^^^-F and in the E'^^^-P plane. 
Moreover, the fraction of bursts detected by GBM which are 
outliers at more than 3cr with respect to the /Speak-^^iso corre- 
lation is ~3%, to be compared with the 6% of outliers found 
(N08) in the BATSE sample. Instead, there are no outliers 
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Table 1. Central values and standard deviations (in brackets) of the distributions of a and E°^^^ for long and short bursts. The table 
also lists the KS probabiUty resulting from the comparison between BATSE and GBM distributions of a and E"^^^. For long bright 
bursts the comparison has been performed by considering (for homogeneity) bursts best modeled by a COMP or a Band model (i.e. 
by excluding from K06 those bursts modeled with a SBPL). 
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Fig. 8. Schematic view of the distribution of long and short GRBs in the E'^^^-F and E'^^^^-P planes. The horizontal dashed line at 
~ 30 keV represents the lower limit for the GBM instrument: the simulations performed in this work show that E'^^^^ can be hardly 
determined below this value. For the BATSE instrument this limit corresponds to ~ 50 keV. The upper limit for BATSE is at ~ 1 
MeV, while for the GBM there is no upper limit in this plane. The vertical dashed line (left panel) shows an example of fluence 
selection, while the dashed curve (right panel) refers to the photon flux selection criterion adopted by G09. 



(at more than 3cr) of the iSpeak - ^iso correlation among GBM 
long GRBs. 

- Short GRBs detected by GBM populate a different region 
in the ^p^^j.-^^ plane with respect to long events, the former 
having larger peak energies and lower fluences compared to 
the latter This is consistent with what found by BATSE and 
confirms that short GRBs do not follow the "Amati" corre- 
lation but they obey the "Yonetoku" correlation defined by 
long events. 

The GBM population of long and short bursts with spec- 
tral information is large enough to allow a statistical compar- 
ison with the BATSE results. For long bursts, we considered 
the fluence distribution of BATSE bursts and we compare it to 
those derived by NIO for GBM bursts. We also compared the 
spectral properties for selected samples of GBM and BATSE 
bursts with well defined E°^\ derived from the spectral analysis. 

peak ^ J 



Two different sample of BATSE bursts are available in litera- 
ture, based on complementary fluence selection criteria. We call 
them faint and bright BATSE samples. We then selected from 
the catalog of NIO two subsamples of GBM bursts based on the 
same fluence criteria applied to the BATSE samples (i.e. 10"^ 
erg/cm^ < F <2x 10"^ erg/cm^ for the faint GBM sample and 
F >2x 10"^ erg/cm^ for the bright GBM sample). 

The £'°''^ distribution derived from the two instruments are 

peak 

quite similar (Fig.[3]and Tab.[T]). Despite its larger energy range, 
the GBM extends the distribution of long bursts only at 
low energies with respect to BATSE. The a distribution, in- 
stead, reveals some difference for the sample of bright bursts: 
BATSE bursts have on average a softer low-energy photon in- 
dex ((acBM) - -0.9 and (ffBATSE) = -1.11, KS probability 
= 6 X lO""*). However, this difference is almost totally due to the 
presence (in the K06 sample of bright BATSE bursts) of GRBs 
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modeled by a smoothly broken power-law (SBPL) function. As 
noted by K06, this model gives a low-energy spectral index sys- 
tematically softer with respect to COMP and Band models. By 
excluding these events, the a distribution of bright BATSE bursts 
is centered around (cfBATSE) = - 1 -00 and the KS probability with 
the GBM is 0.4. 

Also for GBM short bursts we can draw some conclusions 
about their spectral properties. Their distribution is shifted 
towards higher energies compared both to long bursts from the 
same instrument and to short bursts seen by BATSE. The lack of 
low-energy E'^^y^ (below ~ 200 keV) can be accounted for by 
the spectral threshold we derived for the GBM instrument (see 
Fig.|2]i. This hypothesis is supported by the fact that among the 
population of short GBM bursts there are 44 events fitted with 
a curved model (i.e. with E°^\ determined) but there exists a 

^ peak ^ 

large fraction of short bursts (34) whose spectrum is fitted with a 
single power law. On the other hand, the larger energy coverage 
allows the detection of E°^^^, up to ~ 4 MeV. GBM data confirm 
that short bursts have on average a harder a compared to long 
bursts ((aoBM.short) ~ -0.59), as already found in the BATSE 
sample by G09. 

The comparison of short and a representative sample of long 
BATSE GRBs (selected with a similar peak flux threshold) led 
G09 to conclude that their main spectral diversity is due to a 
harder low-energy spectral index in short bursts while their E'^^y^ 
of BATSE is similarly distributed. GBM bursts provide the op- 
portunity of re-examining this result for the population of short 
and long GRBs detected by the GBM and also compare their 
spectral properties with those of the BATSE ones. We find that: 

- E'^^y of short GBM bursts is larger and a smaller that those 
of long one, indicating that short events are harder, both in 
terms of their peak energy and low-energy spectral index. 

- A comparison between GBM and BATSE short bursts re- 
veals that they have similar a while the E'^^^ of short GBM 
bursts is larger than that of short BATSE events (see Fig. |2l 
bottom left panel). This information is allowed by the higher 
energies which can be detected by the GBM. Moreover, the 
diff'erent E°^l^, distribution of BATSE and GBM short bursts 
is affected by the lower sensitivity of the GBM instrument, 
which misses short bursts at low fluences (and therefore low 
E°\). 

- GBM and BATSE long bursts have a similar E'^^^ while 
GBM events tend to have a harder low-energy spectral in- 
dex (Fig.©. 

Fig. [8] shows a schematic representation of the current in- 
formation about the distribution of short and long bursts in the 
E^l^-F and E°^ly-P planes. With respect to BATSE, the GBM 

reveals that long bursts extend to lower E'^^^^, consistently with 
what previously found with other instruments (mainly Hete-II 
and Swift). 

Despite of the high-energy sensitivity, also the ^^^i^ distri- 
bution of GBM long events extends only up to ~ 1 MeV. The 
situation is different for short GRBs whose E'^^^ reach up to ~ 4 

MeV in the present sample. These high £"^^1^ were not detectable 
by BATSE, whose sensitivity drops at ~ 1 MeV (upper horizon- 
tal dashed line in Fig.[8]i. Therefore, the GBM shows that short 
GRBs have larger E'^^^ with respect to long ones, contrary to 
what found with BATSE (G09). 

When comparing the E'^^^ distribution of short and long 
bursts, different conclusions can be drawn, according to the se- 



lection criterion of the samples. The left panel in Fig. [8] shows 
that a given cut in fluence (represented by the vertical dashed 
line) would result in different E°^\ distributions between short 

^ peak 

and long bursts, resulting from their different location in the 
E'^ly-F plane. The right panel in Fig.[8]illustrates, instead, what 
happens for a selection in photon flux. This translates into a curve 
in energy flux: the dashed curve represents the cut applied by 
G09 to select both short and long bursts, corresponding to a 
photon flux larger than 3 ph cm~^s"'. This criterion applied to 
BATSE bursts produces similar E'^^^ distributions of long and 
short events, as indeed found by G09. The very same criterion 
applied to GBM bursts results, instead, in different distributions, 
since short bursts can have very high values, not detected 
in the sample of long bursts. 
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